The Simple Harmonic Oscillator 



Let's investigate the case of no damping (b = 0) a bit more carefully: 

mx" +kx = F ext (t). 

Let's define the parameter 

co n = y/k/m 

and use this to rewrite the equation in terms of m and co n : 

mix" + co 2 n x) = F ext (t). 

We saw this form of the equation in the previous session. Recall that the 
subscript "n" stands for "natural". To remind ourselves why, consider the 
solution in the case of no driving force, i.e. F ext (t) = 0. The characteristic 
equation p(r) = m[r 2 + co 2 ), has roots ±ico n . Thus, the general solution is 

Xh = c\ cos(co n t) + C2 sm(co n t). 

So even without a driver, if we give the system a nudge it will oscillate at 
it's natural frequency oo n . 

Now let's add some sinusoidal input: F ext (t) = Bcos(u>t). Using com- 
plex replacement, we must find a particular solution to 

m(z" + co 2 „z) = Be iu,t . 

Applying the exponential response formula with a = ico, we get 

z = B c iwt = B c iwt . 
p p{ico) m(co1 — co 2 ) 

Taking the real part, 

x„ = 3?(z p ) = —-^ ^cos(a;f) 

y F m{oj 2 — co 1 ) 

is a particular solution. All in all, our general solution is 



x = c\ cos(a;f) + C2 sin(a;f) H -. — , ^- cosicot). 

m[io 2 — co L ) 

Again, the gain is given by 1/ \p(a) \ = l/\p(ico)\. 
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Example. Let's take m = 1, B = 1 and to n = 2, and investigate the resulting 
particular solution as we vary to, the frequency of the driving. The follow- 
ing figure shows the situation for co = 3, 2.5 and 2.1. 




Fig. 1. Solutions for different values of to. 
What do you think happens as to approaches the natural frequency? 

It's no surprise that the solution breaks down when to = to n . This sit- 
uation is called pure resonance and we will investigate it in detail in an 
upcoming session. Notice that it corresponds to the case p(a) = in the 
exponential response formula, since p(ito n ) = 0. For now, we'll note that it 
can be checked that the following is a solution: 

x v (t) = —^—tsm(tot). 

Notice that the extra factor of t before the sine term. So the amplitude grows 
with time, as shown in the following figure. 
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Fig. 2. Solution at pure resonance. 

More generally, the following is the counterpart to the exponential re- 
sponse formula in the pure resonance case, when p(a) =0. 

Resonant Response Formula (RRF). Consider the second order equation 

mx" + he' + bx = Be at , 
with characteristic polynomial p. Then if p(a) = and p'(a) 7^ 0, then 



is a particular solution. 

Once we develop a small amount of algebraic machinery we will be 
able to give a simple proof of this formula. 
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